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ABSTRACT 

A detailed mechanism for the oxidation of benzene Is presented and used 
to compute experimentally obtained concentration profiles and Ignition delay 
times over a wide range of equivalence ratio and temperature. The computed 
results agree qualitatively with all the experimental trends. Quantitative 
agreement Is obtained with several of the composition profiles and for the 
temperature dependence of the Ignition delay times. There are Indications, 
however, that some Important reactions are as yet undiscovered In this 
mechanism. Recent literature expressions have been used for the rate 
coefficients of most Important reactions, except for some Involving phenol. 

The discrepancy between the phenol pyrolysis rate coefficient used In this 
work and a recent literature expression remains to be explained. 

INTRODUCTION 

The major Importance of aromatic hydrocarbons In practical engine fuels 
has made It Imperative to Increase our understanding of the oxidation mechanism 
of these compounds. This knowledge Is necessary for controlling the combustion 
and emission characteristics of all gas turbine combustion systems. The 
simplest aromatic, benzene, has been studied In several recent papers (Venkat 
et al., 1982, McLain et al., 1979, Kern et al . , 1984, Hsu et al . , 1984). A 
review paper on aromatic hydrocarbon oxidation (Brezlnsky, 1986) qualitatively 
outlined a benzene oxidation mechanism and presented a limited number of 
experimentally measured profiles for one benzene-oxygen-nitrogen oxidation in 



a high - temperature turbulent flow reactor. A more comprehensive set of 
concentration profiles for benzene oxidation In the same reactor was recently 
reported by Lovell et al . (1988). 

Ignition delay time measurements (Burcat et al . , 1986) have been reported 
for benzene-oxygen-argon mixtures ignited behind a reflected shock wave. 

These delay times were measured from pressure versus time profiles over a wide 
range of starting conditions. The experiments In the last two mentioned 
papers cover an equivalence ratio range from 0.5 to 2.0 and a temperature 
range from about 1100 to 1600 K. Pressure varied from 1 to about 7 atm. 

This paper presents a comprehensive benzene oxidation mechanism which is 
used to compute the experimental results given In the latter two papers. In 
the following sections we present the results of using the new mechanism to 
compute several concentration profiles for a range of equivalence ratio (<j>) 
from 0.74 to 1.3. Ignition delay times were also computed for equivalence 
ratios of 0.5, 1.0 and 2.0. We show that agreement between experimental and 
computed trends Is generally good. A sensitivity analysis Is also presented to 
demonstrate which reactions have the most significant effect on the computed 
results. All kinetics and sensitivity computations were performed using the 
NASA Kinetics and Sensitivity Code, LSENS (Radhakrl shnan and Blttker, 1986, 
1990). This code Implements the decoupled direct method for sensitivity 
analysis of nonl sothermal systems developed by Radhakrl shnan (1987). 

BENZENE OXIDATION MECHANISM 

The present mechanism is based on the qualitative scheme outlined by 
Brezlnsky (1986), which has been combined with new Information presented in 
several recent papers, as described below. It is also important to use the 
best possible set of hydrogen-oxygen reactions, which are an Important part of 
all hydrocarbon oxidation mechanisms. The hydrogen oxidation scheme developed 
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by Brabbs (1988) was used In these computations. This mechanism was obtained 
by matching experimental Ignition delay times in a shock tube under conditions 
similar to those being used in the present study. A complete listing of the 
reactions and rate coefficients used is given in table I. 

One of the important reactions added to the mechanism of Brezinsky (1986) 
is 

CgHg + C 5 H 5 * C 12 H 10 + H (2) 

It has recently been studied by Fahr and Stein (1988) and had a 
significant impact in obtaining agreement between computed and experimental 
results. As an additional source of H atoms, this reaction increased the 
importance of the main hydrogen-oxygen branching reaction 

H + O 2 * OH + 0 (102) 

Reaction (102) was one of the important steps in the benzene oxidation for all 
the experimental conditions studied in this work. Other significant reactions 
are the attack upon benzene by the radical pool of 0, OH and H species and also 
the reaction of phenyl radical ( CgHs ) with molecular oxygen. These are: 


C 6 H 6 ♦ H * C 6 H 5 + H 2 

(4) 

C 6 H 6 + 0 - C 6 H 5 0 + H 

(5) 

C6 h 6 + 0H * CgHs + H 2 O 

( 6 ) 

c 6 h 5 + o 2 * c 6 h 5 o + 0 

(9) 


As discussed by Brezinsky (1986) and Nicovich et al . (1982), reaction (5) is 
an addition process whose final products could be either those given or else 
phenol. The displacement of a ring hydrogen by oxygen would seem to be a 
simpler procedure than the rearrangement of the adduct to form phenol. 
Therefore, the reaction has been written as the chain-branching process. 

Another important reaction is the decomposition of phenoxy radical 

C 6 H 5 O * C 5 H 5 + CO ( 8 ) 
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The reactions of the cyclopentadlenyl radical, C 5 H 5 , have been written to 
conform to the new experimental results of Lovell et al . (1988). Their results 
for benzene oxidation In the presence of NO 2 Indicated that C 5 H 5 reacts with 
the radical pool (0 and OH) rather than with molecular oxygen. Therefore the 
radical plus C 5 H 5 reactions suggested by Brezlnsky (1986) have been used. 

The formation and destruction reactions of phenol (C 6 H 5 OH) have been 
studied in the recent work of He et al . (1988) and Lovell et al . (1989). The 
latter Investigators studied the phenol pyrolysis In the same flow reactor 
used for the benzene oxidation studies described by Lovell et al . (1988). The 
mechanism and rate coefficients they report were used as the starting point 
for selecting the reactions Involving phenol. This species is formed in the 
benzene oxidation by the reverse of reaction ( 12 ) 

C 6 H 5 OH * C 6 H 5 0 + H (12) 

Reaction (12) exerts an Inhibiting effect on the oxidation and the high rate 
coefficient given by Lovell et al . (1989) strongly supresses the entire 
benzene reaction by competing for H atoms with the H + O 2 chain branching 
reaction. The pre-exponential factor for reaction ( 12 ) had to be reduced to 
one-ninth of Lovell's value. Other reactions of phenol with H and OH were 
taken from Lovell et al . (1989), along with the formation of cyclopentadl ene 


(C 5 H 5 ) by the reaction of C 5 H 5 with phenol. They are 

C 6 H 5 OH + H * C 5 H 5 + OH (13) 

C 6 H 5 OH + H * CgHsO + H 2 (14) 

C 6 H 5 OH + C 5 H 5 * C 6 H 5 0 + C 5 H 6 (15) 

C 5 H 5 OH + OH - CgHsO + H 2 0 (17) 

Phenyl radical dissociates Into two linear molecules according to the reaction 

C 6 H 5 * C 4 H 3 + C 2 H 2 (ID 

A newly measured rate coefficient (Braun-Unkhoff et al . , 1988) for this 


pyrolysis was used here. Their expression Is valid over a wide temperature 
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and pressure range Including the conditions for the data being modeled. The 
measured reaction Is actually the Intramolecular rearrangement of phenyl to 
form "linear" CgHs- This step Is assumed to be followed by the very rapid 
decomposition of that molecule to C 4 H 3 and C 2 H 2 . 

COMPARISON OF COMPUTED AND EXPERIMENTAL RESULTS 
Concentration Profiles 

Lovell et al . (1988) report composition versus time profiles for 
benzene-oxygen-nitrogen oxidations at equivalence ratios of 0.74, 1.0 and 
1.3. The Initial temperature In the flow reactor was 1096 K and remained 
essentially constant. The pressure was 1 atm. The exact zero of reaction 
time Is arbitrary, and was taken as the instant of Injection of the fuel Into 
the hot oxidant stream. In the computations the reactor was modeled as a 
constant pressure (1 atm) homogeneous batch reaction. For all computations In 
this work the thermodynamic data are from the NASA Lewis Research Center data 
base which Is part of the NASA Chemical Equilibrium Composition Computer Code 
of Gordon and McBride (1971). The data for phenyl and phenoxy radicals were 
reported by Burcat, Zeleznlk and McBride (1985). New thermodynamic data for 
several of the minor species were computed by Bonnie J. McBride of NASA Lewis 
Research Center. In particular we should mention that properties for ketyl 
(C 2 HO) and ketene (C 2 H 2 O) were recomputed using a newly computed heat of 
formation for ketyl radical of 38.5 Kcal/mol. New thermodynamic properties 
were also computed for C 5 H 5 , C 5 H 5 O, and C 5 H 4 OH by Dr. A. Burcat. 

Comparisons of computed with experimental profiles are shown in Figs. 1 
and 2. Figure 1 shows benzene and carbon monoxide profiles for the three 
equivalence ratios used and Fig. 2 shows phenol and cyclopentadlene profiles. 
Quantitative agreement between computed and experimental results is generally 
better for benzene and carbon monoxide than for phenol and cyclopentadlene at 
all three equivalence ratios. Figure 1 shows that good quantitative agreement 
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Is obtained for the benzene and carbon monoxide profiles at 4> = 0.74. The 
computed overall rate of the reaction Is slightly slower than the experimental 
one, but the curves have close to the same slope for each species. For 
♦ = 1.0 and 1.3 the computed benzene and CO curves show a somewhat faster 
reaction rate than observed experimentally. The percent difference between 
computed and experimental slopes ranges from 27 to about 50 for these two 
species. This change In the computed net rate from slightly slower to faster 
than the experimental value as equivalence ratio changes is an Indication of 
an Incomplete mechanism. One or more unknown reactions may be unimportant for 
lean mixtures but exert a significant effect on the reaction rate of richer 
mixtures. Possible examples of this could be a reaction of benzene with the 
hydroperoxyl radical or an alternate path for the reaction of benzene with 
molecular oxygen. The phenol concentration profiles in Fig. 2 show that the 
computed curves Increase faster than the experimental ones for all three 
equivalence ratios. The best quantitative agreement Is obtained for <j> = 1.3, 
and all computed curves follow the qualitative trends of the experimental 
lines. The computed curves for cyclopentadl ene show a slower rate of rise 
than do the experimental lines. The final concentration reached by the 
computed curves Is, however, slightly higher than the value for the experimental 
curves. In summary, the computed results match all the experimental results 
qualitatively and there Is some good quantitative matching for 4> = 0.74. 

Ignition Delay Times 

The Ignition delay times, of Burcat et al . (1986) were measured from 
pressure versus time traces obtained by Ignition behind a reflected shock wave. 
The time Interval between shock passage and the first observed "significant" 
pressure rise was taken as the Ignition delay time, x. The reflected shock 
conditions reported by Burcat et al . (1986) were first recomputed with 
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application of a small correction for attenuation of the Initial shock 
velocity. Only thirty-five of the experimental points reported In Burcat 
et al . (1986) were used for our comparisons. All data for equivalence ratio 
of 0.25 were excluded, because examination of the pressure traces showed that 
the pressure rise was very poorly defined for these weak ignitions. In 
addition, all Ignition delay times less than 100 psec were excluded because 
pressure disturbances which cause nonuniform heating of the gas behind the 
reflected shock contribute too high a percentage error to these short ignition 
delay times (Brabbs and Robertson, 1986). To match the experiment, the 
computed ignition delay time, tp, was obtained from each computed pressure 
versus time plot. The computations were performed assuming a constant volume 
batch reaction zone behind the shock wave. For each data point a computed 
Ignition delay time was obtained from the computed pressure profile and was 
defined as the reaction time for a 5 percent rise in pressure. This method 
gave temperature Increases between 30 and 50 K. It was considered to be a 
satisfactory approximation to the experimental technique, which involved the 
reading and interpreting of photographic pressure traces. 

Four mixture conditions among three equivalence ratios were studied. The 
experimental conditions for each mixture are given In Table II. Plots of 
experimental and computed ignition delay time versus the reciprocal of 
temperature for the four mixtures are shown in Figs. 3 to 6. The experimental 
data points are shown for each mixture, as well as a least squares line for 
each set of computed and experimental data. For a given mixture, each set of 
data Is seen to obey an Arrhenius type equation 

t or xp = A exp (AE/RT) 

where R Is the universal gas constant and AE is an apparent activation 
energy which measures the temperature dependence of the delay time for a fixed 
set of Initial concentrations. The agreement between experimental and computed 
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Ignition delay times Is good for mixture 3 ($ = 1.0, strong), fair for mixture 
2 ($ = 1.0, dilute) and poor to fair for the other two mixtures. For $ = 0.5 
and the dilute stoichiometric mixture computed values are seen to be 
consistently higher than the experimental values, whereas the opposite 
situation occurs with the 4> = 2.0 mixture. It Is significant to observe that, 
for all mixtures, the computed Ignition delay times have values of AE quite 
close to the experimental values. A summary of the differences between 
computed and experimental results Is given In Table III. Listed here are all 
the experimental and computed data points and the standard deviations for each 
mixture. In the case of the strong stoichiometric mixture the standard 
deviation of the computed results Is well within the expected experimental 
error. For the lean and the rich mixtures, standard deviations are about 
45 percent, and the direction of the error Is opposite for the two mixtures, 
as noted above. 

In summary, the computations match closely the temperature dependence of 
all the experimental results and quantitatively predict the delay times for 
one stoichiometric mixture. The absolute agreement between computation and 
experiment Is fair to poor for the other mixtures. The most serious 
deficiencies are In predicting delay times of the lean and rich mixtures. For 
the lean mixture the predicted ignition times are too slow, while for the rich 
mixture they are too fast. Similar trends with equivalence ratio were 
obtained for the concentration-profile computations, as mentioned above. These 
results again suggest an unknown reaction. 

Some new measurements of Ignition delay times in ref lected-shock heated 
benzene-oxygen-argon mixtures have recently been reported by Thyagarajan and 
Bhaskaran (1989). Their method of defining ignition delay time Is by detecting 
first light emission, and may or may not be equivalent to the pressure-rise 
criterion of Burcat et al . (1986). These new data agree with those of Burcat 


8 


for some conditions but disagree by a factor of 2 to 3 at other conditions. 

The temperature dependence of the new data Is much weaker than that of the 
Burcat data. The new paper also presents a detailed reaction mechanism which 
computes their experimental Ignition delay time results. However their 
seventy-reaction mechanism contains no steps Involving phenol or cyclopentadlene 
and Is used by the authors to compute concentration versus time profiles only 
for benzene and acetylene. The latter mechanism differs In several ways from 
our mechanism. For example, two Important reactions are written with different 
products. Those given for the benzene plus oxygen atom reaction are phenyl 
and hydroxyl radical rather than phenoxy and H atom used here and recommended 
by Nlcovlch et al . <1982). Also, the products of the phenyl plus molecular 
oxygen reaction are C 2 H 3 , C 2 H 2 and CO, rather than phenoxy and oxygen atom as 
used here and recommended by other Investigators (Brezlnsky, 1986; Lovell 
et al . , 1988). Also no mention Is made in their paper of different colllsional 
efficiencies for the collision partners In the H + O 2 + M recombination 
reaction. For a mechanism to be valid at temperatures around 1000 K, these 
collislonal efficiencies have to be considered. 

In the present work we have used the correlation equation of Thyagarajan 
and Bhaskaran to compute experimental Ignition delay times for several of 
their mixtures and then used our mechanism to compute these same delay times. 
Equivalence ratio was varied from 0.5 to 2 . 0 , pressure from 2 to 10 atm and 
Initial temperature from 1250 to 1750 K. In the worse cases our computed 
results were about 50 percent lower than the experimental values. However, the 
computed delay times agreed within 15 percent or better with the experimental 
results for several cases. Therefore, even though we do not know the exact 
method by which the delay times of Thyagarajan and Bhaskaran are defined, we 
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can say that our present comprehensive mechanism gives computed Ignition delay 
times quite consistent with their experimental data. 

DISCUSSION 

The mechanism presented here successfully explains all the qualitative 
trends and matches quantitatively several concentration profiles and Ignition 
delay times measured during the oxidation of benzene. During the development 
of this mechanism detailed sensitivity analysis computations were performed to 
determine the rate-controlling reactions and to monitor the effect of rate 
coefficient adjustments on the computed species profiles. Figure 7 shows 
normalized sensitivity coefficients for the reactions which control benzene 
concentration In the flow reactor experiments of Lovell et al . (1988). The 
magnitude of a coefficient Is the approximate percent change In the concentration 
which would be caused by a 1 percent change in the rate coefficient of the 
given reaction. A negative sign Indicates that the direction of the 
concentration change Is opposite to the direction of the rate coefficient 
change. Results are shown for the lean and rich equivalence ratio, and a 
reaction time of 50 msec was chosen for this and subsequent figures. Three 
reactions. Independent of equivalence ratio, are seen to be about equally 
Important In controlling the consumption of the fuel; namely reactions (8), 

(12) and (102). The first and last of these reactions promotes the oxidation, 
but reaction (12), which goes In reverse, inhibits the consumption of fuel. 

The Information in Fig. 7 shows the importance of the competing H + O 2 
processes (reactions (102) and (114)), in controlling the oxidation. Both 
reactions accelerate the oxidation process. In the hydrogen-oxygen system, the 
H + O 2 recombination Inhibits the oxidation because of Its chain terminating 
effect. In the benzene-oxygen system, however, hydroperoxyl radical acts as a 
chain carrier by reacting with CgHs radical to form hydroxyl radical (reaction 
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(10)). Figure 7 shows that OH attack on benzene (reaction (6)) Is one of the 
Important steps affecting Its rate of oxidation. Other reactions Involving 
benzene, such as Its pyrolysis and reactions with phenyl and other radicals, 
are also Important In the mechanism. Since we are using a local sensitivity 
analysis method, we can only pinpoint those reactions for which moderate rate 
coefficient changes have an effect on the process. Many other reactions would 
affect the oxidation rate of the fuel If their rate coefficients were very 
uncertain and were changed by large factors. 

Figures 8 to 10 present sensitivity coefficients for carbon monoxide, 
phenol and cyclopentadlene at the same conditions used In Fig. 7. Results are 
similar for all species. The same three reactions control all four 
concentration profiles, and the accelerating effect of the H + O 2 recombination 
can be seen. Reactions (8) and (102) have been studied experimentally and 
their rate coefficients have small uncertainties. We have only adjusted the 
coefficient of reaction (8) by a factor of 1.2 to obtain better CO profiles. 

The rate coefficient of reaction (12), the pyrolysis of phenol. Is not as well 
known. For this process we have used an activation energy equal to the 
endothermlclty of the reaction, 88 Kcal/mol, which agrees with the value used 
by Lovell et al . (1989). The preexponential factor was adjusted to give the 
best possible fit to all the concentration profiles given by Lovell et al . 
(1988). The value obtained, 3.0x10 15 , Is about a factor of nine lower than 
that used by Lovell et al . (1989). This disagreement Is another indication 
that the mechanism may be Incomplete. It should be noted that the Lovell rate 
coefficient Is thermodynamically consistent with the expression given by He 
et al . (1988) for the reverse of reaction (12). However, the lower A factor 
used here for reaction (12) gives similar results when modeling both the flow 
reactor and the shock Ignition experiments. Some experimental data are 
overpredicted and others are underpredicted, depending on equivalence ratio. 
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Sensitivity analysis computations have also been performed to find out 
which reactions were rate-controlling on the Ignition process. Figures 11 to 
13 show sensitivity coefficients of pressure for three different starting 
conditions. All three equivalence ratios are shown and the temperatures used 
range from 1209 to 1435 K. Linder all three conditions, reactions (8), (12) 
and (102) are the most Important ones which affect the Ignition delay time as 
measured by pressure rise. Reactions (8) and (102) accelerate the Ignition 
while reaction (12), which goes in reverse to form phenol, retards the Ignition. 
These are the same reactions that were found to be most sensitive In 
determining the composition profiles for benzene oxidation at lower temperature 
In a flow reactor. These computations show that the H + O 2 recombination 
reaction Is one of the rate controlling steps up to a temperature of 1363 K 
but decreases In Importance as the temperature is raised to 1435 K. As was 
observed for the flow reactor oxidations, this process accelerates the 
Ignition of the fuel . 

This sensitivity analysis shows that It would be difficult to make any 
simple adjustments of rate coefficients In the hope of getting overall better 
agreement with the experimental data. Any changes In rate coefficients would 
Improve agreement at one equivalence ratio, but worsen agreement at a different 
equivalence ratio. One has to search for additional reactions which may resolve 
the discrepancies. These can be reactions among the normal ground-state 
species or possibly reactions among excited-state species. Although the latter 
are not usually considered Important at ordinary combustion temperatures, the 
possibility of excited-state reaction effects cannot be discounted. Reactions 
Involving singlet methylene radical have been suggested as playing a role in 
hydrocarbon oxidation (Miller and Bowman, 1989). However, little accurate 
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Information Is known about the thermodynamics and chemical kinetics of excited 
state species In the benzene-oxygen system, as far as can be determined. So 
the effect of any excited-state species reactions cannot be investigated at 
this time. We have, however, made a series of computations with two reactions 
not used previously. The first one added is the reaction of benzene with 
hydroperoxyl radical 

CgHg + HO2 - C6H5 + H2O2 
AH = 25.6 Kcal/mol 

The activation energy was taken as the heat of reaction and the preexponential 
factor was varied from 1.0x10^ to 1.0x10^. Computations showed that this 
reaction had no effect on any variable profile. In a separate set of 
computations the reaction between benzene and molecular oxygen was written as 
an addition and rearrangement process to give products phenoxy and hydroxyl 
radicals 

C 6 H 6 + 0 2 * C 6 H 5 0 + OH 
AH = 2.3 Kcal /mol 

Both the minimum activation energy of 2.3 Kcal/mol and a higher value of 
10 Kcal/mol were assumed and the preexponential factor was varied over a wide 
range. This reaction had a strong accelerating effect on the oxidation. When 
It was used with the Lovell mechanism for phenol reactions, this reaction 
overcame the strong inhibiting effect of the high phenol reaction rate 
coefficients. However, the oxidation process became too fast, and extremely 
poor concentration profiles and Ignition delay times were obtained. This 
reaction was, therefore, abandoned. 

CONCLUDING REMARKS 

A mechanism for benzene oxidation has been presented which explains 
qualitatively all the observed trends of experimental concentration profiles 
and Ignition delay times over a range of equivalence ratio from 0.50 to 2.0. 

13 


The computations with this mechanism also quantitatively match much of these 
data. However, It Is clear from our comparisons of computed and experimental 
results that this mechanism still must be Improved to obtain agreement between 
computation and experiment and resolve questions about the rate coefficients 
of reactions Involving phenol. The need for additional reactions Is Indicated 
by the fact that the computations underpredict the experimental results for 
some equivalence ratios and overpredict them for other equivalence ratios. 
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TABLE I. - BENZENE OXIDATION MECHANISM 


Number 

Reaction 

A, 

cm 3 , mol , sec 

n 

Ea. 

cal/mol 

Reference and adjustment 
factor (A . F . ) 

1 

CgHg + O 2 * C 5 H 5 + HO 2 

6.30xl0 13 

0 

.0 

60 000 . 

McLain et al . (1979) 

2 

C6 h 6 + c 6 h 5 * Cl2 H 10 + H 

4. OxlO 11 



4 000. 

Fahr and Stein (1988) 

3 

C6 h 6 * C 5 H 5 + H 

5. OxlO 15 



108 000 . 

HSU et al . (1984) (A . F . = 2.0> d 

4 

C 6 H 6 ♦ H = C 6 H 5 + H 2 

3. OxlO 12 



8 100 . 

N 1 cov 1 ch and Ravishankara (1984) 

5 

C 6 H 6 + 0 * C 6 H 5 0 + H 

2.78xl0 13 



4 910. 

Nicovlch et al . (1982) 

6 

C 6 H 6 ♦ OH * C 5 H 5 ♦ h 2 o 

2 . 1 3x 1 O 1 3 



4 580. 

Madronlch and Felder (1985) 

7 

C 4 H 3 + M * C 4 H 2 + H + M 

1.0 xl O 1 6 



60 000 

Miller et al . (1982) 

8 

C 5 H 5 O ’ C 5 H 5 + CO 

2.51x10" 



43 900. 

Lin 

and Lin (1986MA.F. = 1.2) 

9 

C 6 H 5 + 0 2 * C 6 H 5 O + 0 

2 . 1 xl 0 12 



7 470. 

Lin 

and Lin (1987) 

10 

C 5 H 5 + H0 2 * C 5 H 5 O + OH 

5. OxlO 13 



1 000 . 

Estimated 

11 

c 6 h 5 * c 4 h 3 + c 2 h 2 

4.5xl0 13 



72 530. 

Braun-Unkhoff et al . (1988) 

12 

C 6 H 5 OH * C 5 H 5 O ♦ H 

3.0xl0 ]5 



88 000 . 

This work 

13 

C 6 H 5 OH + H * C 6 H 6 + OH 

2 . 2 xl 0 13 



7 910. 

Lovell et al. (1989) 

14 

CsH 5 OH + H * C 6 H 5 0 + H 2 

1 .15xl0 14 



12 400. 

Lovell et al. (1989) 

15 

CfiH^OH + CrHs - CfiH^O + CqHf; 

4 . Ox 1 0 1 4 



25 200. 

This work 

16 

C 5 H 5 + 0 * C 5 H 5 O + H 

1 .OxlO 11 



0 

.0 

Estimated 

17 

C 5 H 5 OH + OH - C 5 H 5 O + H 2 0 

3. OxlO 13 



0.0 

This work 

18 

C 5 H 5 O * C 4 H 5 + CO 

3. OxlO 16 



15 000. 

Estimated 

19 

C 5 H 5 - 0 * C 5 H 5 O 

5. OxlO 13 



0 

.0 



20 

C 5 H 5 + OH * C 5 H 4 OH + H 

1 .OxlO 13 







21 

C 5 H 4 OH * C 4 H 4 + HCO 

1 .OxlO 1 5 







22 

C 5 H 5 + H0 2 * C 5 H 5 O + OH 

2. OxlO 13 






- 

23 

2 C 6 H 5 * Ci 2 H 10 

3 . 1 xlO 1 2 





Colket (1986) 

24 

C 4 H 5 - C 2 H 3 + C 2 H 2 

1 . 4x 1 0 1 3 



32 900. 

Estimated 

25 

c 4 h 2 + 0 * c 2 ho + c 2 h 

1 .OxlO 1 3 



0.0 

McLain (1979) 

26 

c 4 h 2 +• 0 * CO + c 3 h 2 

1 . 2 xl 0 12 



0 

.0 

Miller et al . (1982) 

27 

C 4 H 2 + OH * HCO + C 3 H 2 

3. OxlO 13 



0 

.0 

Miller et al . (1982) 

28 

C 2 H 4 + M * C 2 H 2 + H 2 

9. 33x 10 1 6 



77 200. 

Miller et al . (1982) 

29 

C 2 H 4 + OH # C 2 H 3 + H 2 0 

4.79X10 12 



1 230. 

Westbrook and Dryer (1984) 

30 

C 2 H 4 + 0 * CH 3 + HCO 

3 . 31 xlO 1 2 



1 130. 



31 

c 2 h 4 + 0 * ch 2 o + ch 2 

2 . 51 x 10 1 3 



5 000. 



32 

C 2 H 4 + OH * CH 3 + CH 2 0 

2. OxlO 12 



960. 



33 

C 2 H 3 + M * C 2 H 2 + H + M 

7 . 94x1 O 1 4 



31 500. 

Westbrook and Dryer (1984XA.F. - 0.38) 

34 

C 2 H 3 + 0 2 - CH 2 0 + HCO 

3.98xl0 12 



-250. 

Slagle and Gutman (1986) 

35 

C 2 H 3 + H * C 2 H 2 + H 2 

6 . OxlO 12 



0 

0 

Miller et al . (1982) 

36 

C 2 H 3 + OH * C 2 H 2 + H 2 0 

5. OxlO 12 







37 

C 2 H 3 t ch 2 * c 2 h 2 + ch 3 

3-OOxlO 1 3 







38 

C 2 H 3 * C 2 H * 2 C 2 H 2 

3 .OOxlO 1 3 







39 

C 2 H 3 + 0 * C 2 H 2 0 + H 

3.3x10 13 







40 

ch 2 ♦ ch 2 * C 2 H 2 -*■ H 2 

4. OxlO 13 







41 

CH 2 + CH 2 * C 2 H 3 * H 

5. OxlO 12 





Westbrook and Dryer (1984) 

42 

CH 2 + OH * CH ♦ H 2 0 

2 . 5 1 x 1 0 1 1 

0 . 

67 

25 700. 

Miller et al . (1982) 

43 

CH 2 + 0 * CH ♦ OH 

2 . 0 x 10 ' 1 

0 . 

68 

25 000. 



44 

CH 2 + 0 2 * C0 2 + 2H 

1 . 59x 1 O 1 2 

0 . 

0 

1 000 . 



45 

C 2 H 2 ♦ M * C 2 H 4 H 

4. 1 7x 1 O 1 6 



107 000. 



46 

C 2 H 2 + C 2 H 2 * C 4 H 3 + H 

2. OxlO 12 



45 900. 



47 

C 2 H 2 + 0 * CH 2 + CO 

1 . 6 x 1 0 1 4 



9 890. 

Frank et al . (1 986) 

48 

C 2 H 2 + 0 * C 2 HO + H 

4. OxlO 14 



10 660. 

Frank et al . (1986) 

49 

C 2 H 2 + OH * C 2 H + H 2 0 

6. 31 xlO 12 



7 000 

Westbrook and Dryer (1984) 

50 

C 2 H 2 + OH s C 2 H 2 0 + H 

3.2x10 ' 1 



200 . 

Westbrook and Dryer (1984) 

51 

C 2 H 2 + C 2 H * C 4 H 2 + H 

3. OxlO 13 



0 . 

0 

Miller et al . (1982) 

52 

C 2 H 2 + CH 2 - C 3 H 3 + H 

1 . 2 x 10 13 



6 600. 

Bohland et al . (1 986) 

53 

C 3 H 4 + M - C 3 H 3 + H + M 

2. OxlO 17 



65 000. 

Pamldlmukkala et al . (1987) 

54 

C 2 H 2 0 + OH * CH 2 0 f HCO 

2 . 8 xl 0 13 



0 . 

0 

Miller et al . (1982) 

55 

C 2 H 2 0 + OH * C 2 HO + H 2 0 

7.5x10 12 



3 000. 

Miller et al . (1982) 

56 

C 2 H 2 0 + H * CH 3 + CO 

1.13xl0 13 



3 428. 

Westbrook and Dryer (1984) 

57 

C 2 H 2 0 4- H * C 2 HO + H 2 

7.5xl0 ' 3 



8 000 . 

Miller et al. (1982) 

58 

C 2 H 2 0 + 0 * C 2 HO + OH 

5. OxlO 13 



8 000 . 

j 


59 

C 2 H 2 0 + 0 * CH 2 0 + CO 

2. OxlO 13 



0 . 

Li 

i 
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- Concluded. 


Number 

Reaction 

60 

C 2 H 2 0 + M * CH 2 + CO + M 

61 

C 2 HO + 0 2 - 2 CO + OH 

62 

C 2 H0 + 0 * 2 CO + H 

63 

C 2 H0 + OH * 2 HCO 

64 

C 2 H0 + H * CH 2 4 CO 

65 

C 2 H0 + CH 2 * C 2 H 3 + CO 

66 

C 2 HO + CH 2 * CH 2 0 + C 2 H 

67 

c 2 ho + c 2 ho * C 2 H 2 «. 2 CO 

68 

C 2 H 4 OH * C 2 HO 4 H 

69 

c 2 h + 0 2 * C 2 H0 4 0 

70 

C 2 H + 0 - CO + CH 

71 

CH 4 4 M a? CH 3 4 H 4 M 

72 

CH 4 + 0 2 * CH 3 4 ho 2 

73 

CH 4 4 H ^ CH 3 4 H 2 

74 

CH 4 4 OH ^ CH 3 4 h 2 o 

75 

CH 4 4 0 ^ CH 3 4 OH 

76 

CH 3 + 0 2 * CH 3 0 4 0 

77 

CH 3 4 OH 5 CH 3 0 4 H 

78 

CH 3 0 + M * CH 2 0 4 H 

79 

CH 3 0 4 0 2 * CH 2 0 4 ho 2 

80 

CH 3 0 4 H * CH 2 0 4 H 2 

81 

CH 3 4 CH 3 * C 2 H 4 4 h 2 

82 

CH 3 4 0 * CH 2 0 4 H 

83 

CH 3 4 CH 2 0 * CH 4 4 HCO 

84 

CH 3 4 HCO * CH 4 4 CO 

35 

CH 3 4 H0 2 * CH 3 0 4 OH 

86 

CH 2 0 4 M - HCO 4 H 4 M 

87 

CH 2 0 4 OH ^ HCO 4 H 2 0 

88 

CH 2 0 4 H - HCO 4 H 2 

89 

CH 2 0 4 0 ^ HCO 4 OH 

90 

HCO 4 H0 2 * CH 2 0 4 0 2 

91 

HCO 4 M * H 4 CO 4 M 

92 

HCO 4 0 2 * CO 4 H0 2 

93 

HCO 4 OH * CO 4 H 2 0 

94 

HCO 4 H * CO 4 H 2 

95 

HCO 4 0 ^ CO 4 OH 

96 

CH 4 0 2 * HCO 4 0 

97 

CO 4 0 4 M * C0 2 4 M 

98 

CO 4 0 2 - C0 2 4 0 

99 

CO 4 OH * C0 2 4 H 

100 

CO 4 H0 2 * C0 2 4 OH 

101 

0 4 H 2 0 * OH 4 OH 

102 

H 4 0 2 * OH 4 0 

103 

0 4 H 2 * OH 4 H 

104 

H + H0 2 * H 2 4 0 2 

105 

0 4 H0 2 ^ OH 4 0 2 

106 

H0 2 4 OH % H 2 0 4 0 2 

107 

H 4 H0 2 - OH 4 OH 

108 

H 2 4 H0 2 * H 2 0 2 4 H 

109 

OH 4 H 2 0 2 * H 2 0 4 H0 2 

110 

H0 2 4 H0 2 # H 2 0 2 4 0 2 

111 

H 4 H 2 0 2 * OH 4 H 2 0 

112 

H 2 0 2 4 M * OH 4 OH 4 M 

113 

H 2 4 OH * H 2 0 4 H 

114 

H 4 0 2 4 M * H0 2 4 M 

115 

H 2 0 4 M * H 4 OH 4 M 

116 

H 4 0 4 M * OH 4 M 

117 

H 2 4M*H4H4M 

118 

0 2 4 H ’ 0 4 0 4 M 


A, 

cm3, moi ? sec 

2 . Oxl 0 1 6 
1 ,46xl0 12 
1 . 20 xl 0 12 
1 .OxlO 13 
5.0x10 13 
3. OxlO 13 
1 .OxlO 13 
1 .OxlO 13 
2. OxlO 13 
5. OxlO 13 
5. OxlO 13 
2. OxlO 17 
7 . 94x1 O 1 3 
1 .26xl0 14 
2 . SOxlO 1 3 
1 .90xl0 14 
4 . 7 9 x 1 0 1 3 
6 . 30x 1 0 1 2 
5. OxlO 13 
1 .OxlO 12 
2 . OxlO 13 
1 .OxlO 16 
1 . 2 9x 1 0 1 4 
1 .OxlO 10 
3. OxlO 11 
2 . 00 x 1 0 1 3 
5. OxlO 16 
3. OxlO 13 
2 . 5x JO 1 3 
3.5xl0 13 
1 .OxlO 14 
2.94X10 14 
3 . 31 xlO 12 
1 .OxlO 14 
2. OxlO 14 
1 .OxlO 14 
1 .OxlO 13 
5 . 9x 1 0 1 6 
2.5xl0 12 
4 . HxlO 1 1 
5.75x10 13 
6 . 8 x 1 O 1 3 
1 . 89x1 0 1 4 
4 . 20x 1 0 1 4 
7 . 2 8 x 1 0 1 3 
5. OxlO 13 
8, OxlO 12 
1 .34xl0 14 
7 . 9 1 x 1 0 1 3 

6. 1 xlO 1 2 
1 .8xl0 12 
7 . 8 x 1 0 1 1 

1 . 44x1 O 1 7 
4 . 74x1 0 1 3 
1 . 46x1 0 1 5 
1 .30xl0 15 
7 . 1 x 1 0 1 8 

2 . 2 x 1 0 1 4 
1 .8xl0 18 


Reference and adjustment 
factor (A . F . ) 


60 000 . I 

2 500. 1 

0.0 Westbrook and Dryer (1984) 
t Ml 1 ler et a! . (1982) 


2 000 . 
0.0 
0.0 

1 500. 

0.0 

88 000 . 
56 000. 

11 900. 

5 010. 

11 720. 
29 000. 

0.0 

21 000 . 

6 000 . 

0.0 

32 000. 

2 000 . 

6 000 . 

0.0 

0.0 

76 500. 

1 200 . 

3 990. 

3 510. 

3 000. 
15 570. 

7 000. 
0.0 


Westbrook and Dryer (1984) 
Westbrook and Dryer (1984) 
Brabbs and Brokaw (1974) 
Brabbs and Brokaw (1974) 
Brabbs and Brokaw (1974) 
Westbrook and Dryer (1984) 
Westley (1980) 

Westbrook and Dryer (1984) 


Warnatz (1984) 


4 100. 
47 690. 

1 000. 
22 930. 
18 365. 
16 400. 
13 750. 

2 126. 

1 000 . 
0.0 
1 070. 
25 000. 

1 430. 
0.0 
0.0 

45 510. 

6 098. 
-1 000. 
105 140. 
0.0 

96 000. 
118 000 . 


Westley (1980) 

Cherian et al . , (1981)° 
Westbrook and Dryer (1984) 


Brabbs and Brokaw (1974) 
Westbrook and Dryer (1984) 
Brabbs and Muslak (1988) 


Adjustment for Increased pressure (approximate) from Kiefer et al . (1985) used for Ignition-delay computations 
^Computed from reverse reaction rate coefficient and equilibrium constant. 

COLLISIONAl EFFICIENCIES 


Reaction 112 
Reaction 1 14 
Reaction 115 
Reaction 117 


H? = 2.3; 0? = 0.78; H?0 - 6.0; H 2 0 2 - 6-6 

hb = 3.0; 0 2 - 1.3; H 2 0 - 21.3; N 2 - 1.3; C0 2 - 7.0; C 6 H 6 - 20.0; CH 4 = 5.0 

H 2 » 4.0; 0 2 - 1.5; H 2 0 « 20.0; N 2 = 1.5; C0 2 * 4.0; C 6 H 6 * 20.0 

H 2 - 4.1; 0 2 - 2.0; H 2 0 * 15.0; N 2 - 2.0 


[Rate constant used equals tabulated expression multiplied by adjustment factor.] 




TABLE II. - INITIAL MIXTURE CONDITIONS FOR IGNITION DELAY TIME MEASUREMENTS 


Mi xture 

Equivalence 

ratio, 

♦ 

Benzene , 
mol t 

Oxygen, 
mol % 

Argon, 
mol 1 

Initial 

temperature 

range, 

K 

Initial 
pressure , 
atm 

i 

0.5 

1.354 

20.313 

78.333 

1209-1345 

1 .9-2.2 

2 

1 .0 

.516 

3.868 

95.616 

1345-1528 

5.6-7. 1 

3 

1 .0 

1.690 

12.675 

85.635 

1283-1435 

2. 0-2. 5 

4 

2.0 

1.354 

5.093 

93.553 

1363-1600 

2. 0-2. 6 


TABLE III - COMPARISON OF COMPUTED AND EXPERIMENTAL INGINTION DELAY TIMES 


Mixture number 
and description 

Initial 
temperature , 
K 

Experimental 
ignition 
delay time, 
psec 

Computed Ignition 
delay time, 
psec 

Percent 

difference 

Percent 

standard 

deviation 

1. Equivalence 

1209. 

878. 

1200. 

36.7 

44.5 

ratio = 0.5 

1 227. 

743. 

960. 

29.2 

— 


1254. 

435. 

600. 

37.9 

— 

p = 2 ATM. 

1 276. 

330. 

480. 

45.4 

— 


1291 . 

272. 

390. 

43.4 

— 


1307. 

185. 

315. 

70.3 

— 


1314. 

202. 

300. 

48.5 

— 


1 345. 

159. 

209. 

31.4 

— 

2. Equivalence 

1345. 

755. 

640. 

-15.2 

24.5 

ratio =1.0 

1 374. 

604. 

445. 

-26.3 

— 

( d 1 lute) 

1402. 

415. 

332. 

-20.0 




1412. 

412. 

295. 

-28.4 

— 

p £ 5-7 ATM. 

1 428. 

367. 

250. 

-31.9 

— 


1482. 

213. 

147. 

-31.0 

— 


1 525. 

122. 

100. 

-18.0 

— 


1 528. 

122. 

99. 

-18.9 

— 

3. Equivalence 

1283. 

750. 

760. 

1.3 

8.9 

ratio =1.0 

1290. 

613. 

700. 

14.2 

— 

(strong) 

1294. 

607. 

660. 

8.7 

— 


1 328. 

490. 

440. 

-10.2 

— 

p s 2 ATM. 

1 355. 

303. 

320. 

5.6 

1 


1369. 

287. 

280. 

-2.4 

— 


1 379. 

291 . 

258. 

-11.3 

— 


1405. 

198. 

208. 

5.1 

— 


1408. 

189. 

178. 

-5.8 

— 


1417. 

178. 

163. 

-8.4 

— 


1435. 

151 . 

130. 

-13.9 

— 

4. Equivalence 

1363. 

1520. 

780. 

-48.7 

45.6 

ratio = 2.0 

1415. 

890. 

440. 

-50.6 

— 


1457. 

599. 

320. 

-46.6 

— 

p = 2 ATM. 

1540. 

274. 

150. 

-45.3 

— 


1554. 

243. 

130. 

-46.5 

— 


1570. 

211 . 

110. 

-47.9 

— 


1582. 

157. 

100. 

-36.3 

— 


1600. 

154. 

90. 

-41.5 

— 
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CONCENTRATION, PPM (VOLUME) CONCENTRATION, PPMV CONCENTRATION, PPMV 


EXPERIMENTAL 

COMPUTED 


EXPERIMENTAL 
COMPUTED > 



(a) <j> m 0.74. 



(b) q>= 1.0. 



0 20 40 60 80 100 1 20 140 1 60 1 80 

TIME, msec 
(c) tp = 1 .3. 

Figure 1 . - Benzene and CO concentration 
versus time, T 0 = 1096 K. 


z _ 
O UJ 

160 


120 

S3 

80 

02 
Z CL 

40 

O Q- 


O 

0 


C 6 H 5 OH 



20 40 60 80 100 120 140 160 180 

TIME, msec 
(c) q> = 1.3. 

Figure 2. - Cyciopentadiene and phenol 
concentration versus time, T 0 = 1096 K 


□— EXPERIMENT 

COMPUTED 


§> 1000 

3 . 



7.2 7.4 7.6 7.8 8.0 8.2 

i n rton/T ^"1 


1350 1280 1219 

TEMPERATURE, K 

Figure 3. - Benzene-oxygen-argon ignition 
delay time versus reciprocal of 
temperature, <p = 0.5. 



IGNITION DELAY TIME, nsec 



10 000/T, K _1 


1515 1430 1350 

TEMPERATURE, K 

Figure 4. - Benzene-oxygen- argon ignition 
delay time versus reciprocal of temper- 
ature, (p - 1 .0, dilute mixture 95.6% Ar. 



10 000/T, K" 1 


1430 1350 1280 

TEMPERATURE, K 


Figure 5. - Benzene-oxygen -argon ignition 
delay time versus reciprocal of temper- 
ature, 9 » 1 .0, strong mixture 85.6% Ar. 



1615 1515 1430 1350 

TEMPERATURE, K 

Figure 6. - Benzene-oxygen-argon ignition 
delay time versus reciprocal of 
temperature, <p = 2.0. 
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(a) Equivalence ratio 0.74, temperature 1096 K, time = 50 msec. 



(b) Equivalence ratio 1 .3, temperature 1096 K, time = 50 msec. 


Figure 7. - Sensitivity coefficients for benzene concentration. 
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No. 


REACTION 


SENSITIVITY COEFFICIENT 


-1.4 -1.2 —1.0 -0.8 -0.6 -0.4 -0.2 0 .2 .4 .6 .8 1.0 1.2 1.4 


I I I I I I I I 


8 c 6 h 5 o^c 5 h 5 + co 

12 C 6 H 5 OH ^ CgHgO ♦ H 

102 H + 0 2 ^0H + 0 

114 H + 0 2 + M^H0 2 +M 

5 CgHg + O ^ C^O ♦ H 

17 C5H5OH + OH ^ CgH 5 0 ♦ H2O 
1 9 C^H 5 + O CgHgO 



(a) Equivalence ratio 0.74, temperature 1096 K, time * 50 msec. 


No. 

REACTION 

SENSITIVITY COEFFICIENT 

-1.4 -1.2 -1.0 -0.8 -0.6 -0.4 -0.2 
1 1 1 1 1 1 1 


8 

C 6 H 5 0 — C 5 H 5 + CO 


S35SS3333531H 

12 

c 6 h 5 oh^c 6 h 5 o + h 



102 

h + o 2 ^oh + o 


5 

114 



1 


CgHg + O^CeHsO + H 







19 

CsHg + O^CgHsO 


9HI 


(b) Equivalence ratio 1 .3, temperature 1096 K, time = 50 msec. 
Figure 8. - Sensitivity coefficients for carbon monoxide concentration. 
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REACTION 


SENSITIVITY COEFFICIENT 


No. I 


17 


CgHgOH + OH ^ C e H s O + 


12 I 
8 

1021 

6 

20 | 


C 6 H 5 OH ^ C 6 H 5° + H 


CgHgO^CsHs^CO 


h + o 2 ^oh + o 


CfiHe + OH ^ CgHs + h 2 o 


C5H5 + OH ^ C5H4OH + H 


114 


H + 0 2 + M ^ H0 2 + M 



(a) Equivalence ratio 0.74, temperature 1096 K, time = 50 msec. 


NO. 


SENSITIVITY COEFFICIENT 

REACTION 

- 8 -.6 -.4 -.2 C 

1 1 1 1 1 

.2 .4 .6 .8 

1 1 1 1 

102 

h + o 2 ^oh + o 


1 

12 

C'6 H 5 OH ^ C 6 H 5° * H 

5 

— 

8 

C 6 H 5° ^ C 5 H 5 + CO 



D 


xz/y/zzi 


114 



s 

20 

C 5 H 5 + OH ^ C 5 H 4 OH + H 


| 

6 

CeHe + oh ^ CgHs + h 2 o 




1 

(b) Equivalence ratio 1 .3, temperature 1096 K, time = 

50 msec. 


Figure 9. - Sensitivity coefficients for phenol concentration. 
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No. 


REACTION 


SENSITIVITY COEFFICIENT 


102 H + Oj^OH + O 

8 C 6 H 5 0 it CgHg + CO 

12 CgHgOH it CgHgO + H 

15 CgHgOH + CgHg CgHgO + CgHg 
17 CgHgOH + OH it CgHgO + HjO 
® CgHg + OH CgHg + H 2 O 
19 CgHg + O it CgHgO 



(a) Equivalence ralio 0.74, temperature 1096 K, time - 50 msec. 
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Figure 12. - Sensitivity coefficients for pressure; equivalence ratio 1.0 (85% Ar), temperature 1435 K, 
time = 100 psec. 
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No. 


REACTION 


SENSITIVITY COEFFICIENT 


-.04 -.03 -.02 -.01 0 .01 .02 .03 04 

J I I I | I I I 1 


102 


12 


h + o 2 ^oh + o 


c 6 h 5 oh^c 6 h 5 o+ h 


xzzzzzzzzzzzzz 




C 6 H 5 0 ^ C5H5 + CO 


C 6 H 6 + O ^ CgH^O + H 


WWWWWM 


fttz 


19 

[114 


C5H5 + O ^ C5H 5 0 


SS3 


H + 0 2 + M^H0 2 + M 


K3 


17 


CgH 5 OH + OH ^ CgHgO ♦ H 2 0 




Figure 13.- Sensitivity coefficients for pressure; equivalence ratio 2.0, temperature 1 363 K. time = 600 nsec. 
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